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Abstract. We have re-analyzed using a hydrodynamic model large flaring events on three different 
categories of pre-main sequence (PMS) stars: the young stellar object (YSO) YLW 15, the classical 
T Tauri star (CTTS) LkHQ92, the weak-line T Tauri star (WTTS) V773 Tau, and the WTTS 
HD 283572 (the first three objects were observed by ASCA, the last by ROSAT; all have been 
previously reported in the literature). The first three flares were previously analyzed on the basis 
of the quasi-static model mostly used up to now, consistently yielding large loops (L i?,) and no 
evidence of sustained heating. Our hydrodynamic modeling approach, however, shows that the size 
of the flaring regions must be much smaller (i ;^ i?*) and moreover this method shows in all cases 
evidence of vigorous sustained heating during the flare decay, so that the decay of the observed 
light curve actually reflects the temporal profile of the heating rather than that of the free decay of 
the heated loop(s). The events on the protostar YLW 15 have durations comparable to the stellar 
rotation period, so that their limited size and their lack of self-eclipses give evidence of a polar 
location on the star. This is in contrast with the recently advanced hypothesis that these flares 
are due to long loops spanning the region between the star and the accretion disk. In general, the 
present analysis shows that flaring coronae on PMS stars have a structure similar to the coronae 
on older active stars. 
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1. Introduction 

Strong X-ray emission is associated with all stages of the 
early evolution of low-mass stars, from the protostellar 
Young Stellar Objects - YSO or Class I - to the Classical 
T Tau - CTTS or Class II - to the Weak-Line T Tau phase 



WT"i S or Class ill. i he salient characteristics ot stel- 
lar activity m the pre-mam sequence (^FiVlS) phase have 



been recently reviewed by Feigelson fc Montmerle (1999| ) : 
in most cases the observed X-ray emission has been inter- 
preted as due to magnetically confined (and likely mag- 
netically heated) plasma, although with much enhanced 



and CTTS's significant amounts of circumstellar material 
(mostly in the form of a disk) is present, and magnetic 
fields are thought to fun nel and modulate the accretion 
onto the (proto)star (e.g. [llartmann, 1998 ). Various pos- 
sibilities for the magnetic interaction between th e star and 
the disk are for example shown in Figs. 2 and 3 of Feigelson 
fc Montmerle (1999| ); understanding the topology and di- 
mensional scales of the coronal plasma (and, indirectly, of 
the magnetic field) can, in the case of PMS stars, give im- 
portant constraints to the accretion process. The stellar 
activity also has a significant influence on the circumstel 
lar environment (see e.g. the review of Glassgold et al. 



with respect to mam-sequence stars of comparable mass 
Whether the corona is however just an enhanced version of 
a "solar-like" one, or whether significantly different mech- 
anisms (and coronal topologies) are at work has been a 
matter of debate, in particular in relationship to the possi- 
ble influence of the accretion disk on the corona: in YSO's 



both | m absolute and relative terms ^ activity levels 2OOOD , and for example a low-lying X-ray emitting region 
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will irradiate (and thus ionize) the circumstellar disk much 
less than a very extended emitting region from which X- 
rays can reach the disk out to much larger distances from 
the parent star. 

The study of the decay of flares has in the last decade 
been extensively used to model the characteristics of the 
flaring regions in different stellar types, and in particular 
their size. This in turn has been used to derive constraints 
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on the structuring of the active corona, with results which 
are however, as discussed below, model-dependent. The 
quasi-static method of flare modeling (originally described 
by van den Oord & Mewc, f 989) assumes that the flar- 
ing loops decay freely after an initial impulsive 
heating event, and that the decay can be described 
through a series of static states; while in principle the 
formalism allows for the presence of heating during the de- 
cay phase, in practice in no case has the application of the 
formalism resulted in sustained heating being detected. 
Thus, in the presence of the long decay times typically 
observed for large flares in active stars, the quasi-static 
method consistently yields long loops and extended coro- 
nae. This method has been extensively applied in the last 
decade, resulting in a general assumed framework of coro- 
nae of active stars populated with very long, extended loop 
structures. In the case of PMS stars, these have often been 
thought to extend to, and link with, the accretion disk, re- 
sulting in a distinctly "non-solar" coronal geometry. 



More recently, Reale et al. (1997 ) have developed a 
method based on hydrodynamic modeling of decaying flar- 
ing loops, which is able to detect the presence of heating 
during the decay and to adjust the estimate of the loop 
length accordingly, and which has been tested ("cali- 
brated") with good results on solar flares of which 
the loop length can be checked from images. When 
applied to the only stellar flare for which an eclipse al- 
lows to make a geometric size estimate (the Algol SAX 



flare, [Schmitt fc Favata, 1999| ; [Favata fc Schmitt, 1999| ) 



the quasi-static method has been shown to over-estimate 
the loop size by a large factor, while the hydrodynamic 
modeling diagnoses the presence of long-lasting sustained 
heating during the decay, and yields a more realistic esti- 
mate of the region's size (although still overestimating it). 
Application of this method to a number of events on dif- 
ferent types of coronal sources has shown that sustained 
heating is an essential feature of large stellar flares, and 
that therefore the flaring structures are actually signifi- 
cantly smaller than so far postulated. This applies to flare 
stars ( CM L eo, [Reale fc Micela, 199^; E V Lac, [Favata 
et al., ^OOOct AD Leo, [Favata et al., 2000a|), to young sin- 
gle solar- type stars (AB Dor, Maggio et al., 20()0[ ), to ac- 
tive bi naries (AR Lac and CF Tuc, Favata, 2000) an d to 
Algol ([Favata fc Schmitt, 19991 ; [Favata et al., 2000b[) . 

The hydrodynamic modeling approach is based on de- 
tailed simulations of decaying flaring loops. The simula- 
tions show that the key diagnostic is the slope C of the 
flare's decay in the logT vs. logn plot (actually, in the 
case of stellar flares, where no direct d ensity diagnostics 
are usually available, the logT vs. log -\/ EM^ plot): if sus- 
tained heating is present, the slope will be shallower than 
in the case of an impulsively heated, freely decaying loop 
(cf. e.g., Fig. |l|). A fit to the parameters of the numeri- 
cal model allows to derive the length of the flaring loop as 
a function of the observed decay time scale tlc , of the ( 



parameter and of the peak temperature in the flaring loop 
Tmax- The latter quantity is defined as the maximum 
intrinsic temperature of a model loop whose spectrum syn- 
thesized in the band of a given X-ray spectrometer yields, 
when fit with a single-temperature spectrum, a maximum 
observed temperature Tp. 

Large flares on PMS stars have in the past always been 
studied with the quasi-static method, resulting in typical 
sizes of L ~ 10^^-10^^ cm, i.e. L > i?*, which have often 
lead to the hypothesis of loops linking the photosphere 
with the accretion disk (e.g. Tsuboi et al., 2000| ). In the 
present paper we study, through the hydrodynamic mod- 



eling method of Reale et al. (1997), stellar X-ray flares 
observed on PMS stars from Class I to Class III, with the 
aim of determining whether they feature strong sustained 
heating (as observed in all stellar types studied so far with 
this approach) and thus whether the related coronal struc- 
tures are significantly more compact than postulated up 
to now. Given the relatively large average distance of star- 
forming regions (and therefore of PMS stars) from Earth, 
only few (and very intense) X-ray fiares from PMS stars 
have been observed with sufficient statistics, and therefore 
the choice of the sample studied in this paper is driven by 
the (random) availability of the data, and is neither com- 
plete or unbiased. 

The present paper is so structured: Sect. || gives a sum- 
mary description of the sample's characteristics, in Sect. |^ 
the analysis of each flare is discussed in detail, while in 
Sect. ^ the results are discussed and summarized. 



2. Sample studied 

We have identified, in the literature, four X-ray fiares from 
PMS's which have sufficient statistics to allow their study 
with the hydrodynamic modeling approach: a flare ob- 



served with ASCA on the protostar YLW 15 (Tsuboi et al 
^000[ ), a flare observed with the ROSAT PSPC on the 



CTTS LkHa 92 ( [Preibisch et al., 1993|), a fiare observed 
with ASCA on the WTTS V773 Tau ( [Tsuboi et al., l"998|) 
and a PSPC fiare on the WTTS HD 283572 ( [Stelzcr et al 



^000[) . The sample (although small) spans the complete 



pre-main sequence phase and covers all key evolutionary 
stages; the somewhat older star AB Dor, just arriving on 
the main sequence, has been studied with the same ap- 
proach by Maggio et al. (2000 ), thus linking the present 
sample with the subsequent evolutionary stage. 

The first three events discussed here have been previ- 
ously analyzed with the quasi-static method, while the 



fourth has been reported by Stelzer et al. (200C) but 
not analyzed in detail. The quasi-static analysis has re- 
sulted (as typical for large stellar flares) in large loop 
sizes (L ~ 3.5 i?* for YLW 15, L ~ 0.6i?* for LkHa 92, 
i ~ 1.3 i?, for V773 Tau). 



^ where the emission measure is defined as EM = 

J rieriudV 



3. Analysis of the flaring events 
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Fig. 1. The evolution of the YLW 15 flare. The top panel 
shows the evolution of the flare's emission measure (used 
as a proxy to the light curve, e.g. to the X-ray count rate), 
together with the best-fit exponential decay, while the bot- 
tom panel shows the evolution of the flare's decay in the 
logT vs. logn plane. The solid line is the best fit to the 
decay {( = 1.0), while the dotted hue is the path the de- 
cay would have followed if it had been quasi-static with 
no sustained heating during the decay phase (Cgs — 1-7). 



3.1. YLW 15 

YLW 15 is a low-mass Class I protostar in the p Oph star- 
forming region (d ~ 165 pc), w ith an estimated radiu s 
i?, ~ 4.2 Rq = 2.9 X 10" cm (e.g. [Montmerle et al., 20001) . 
Its high X-ray activity was noti ced through a large flare 
observed with the ROSAT HRI ( |Grosso ct al., 199?!) ; the 
lack of spectroscopic capabilities of the HRI however pre- 
vents a detailed analysis of this event. YLW 15 was later 
observed with ASCA, which detected three consecutive 
significant flaring events, which are discussed and ana- 
lyzed by Tsuboi et al. (2000 ). The first and most intense 
of the three events has a decay time-scale tlc = 31 ks, a 
measured maximum temperature Td ~ 64 MK and 
a peak X-ray luminosity of ~ 2 x lO'^^ 



erg s 



Tsuboi 



et al. t|2000| ) have analyzed the event using the quasi 
static framework (assuming impulsive heating), deriving 
a semi-length of the flaring loop L ~ 10^^ cm ~ 3.5 i?*, 
a density n ~ 5 x 10^" cm~^ and a minimum value, for the 
magnetic field, of -B ~ 150 G. Similar physical parameters 
are derived for the second and third event (which have 
however lower statistics). 

Using the spect ral parameters {T and EAI) derived by 
Tsuboi et al. (200C ) for three temporal intervals during the 
flare decay we have modeled the first flare using the hydro- 
dynamic modeling approach of Reale et al. (1997 ), recal- 
ibrated for use with the ASCA detectors as d escribed by 
Favata et al. (2000c] ) and [Favata et al. (2000a| ). Although 



the flare decay has been divided in only three time in- 
tervals (due to the limited number of photons detected), 
this is sufficient to constrain the presence of sustained 
heating and thus the loop's length. The flare evolution 
and its best-fit decay in the logT vs. logn plane are 
shown in Fig. 0. The slope of the decay, C = 1.0 ± 0.3 
is significantly shallower than the slope that would be ob- 
served if the decay were indeed quasi-static (Cgs = 1.7). 
The corresponding ratio between the the observed de- 
cay time and the intrinsic thermodynamic decay cooling 



time of the loop without additional heating ( 3erio et al _ 
199l|) is TLc/Tth = F(C) ~ 2.7[2.2 - 3.6]|, indicating the 



presence of significant sustained heating during the de- 
cay. The actual peak temperature of the plasma in the 
flaring loop is Tmax — 150 MK, and the loop length is 
L ~ 4.6[3.5 - 5.6] x 10"cm ~ 1.6[1.2 - 1.9] i?*, a factor 
of ~ 2 smaller than estimated by the quasi-static method 
under the assumption of free decay. 

To estimate the density n of the plasma and the 
strength of the confining magnetic field B, the volume 
V of the flaring region and therefore its geometry must 
be known. Assuming a loop with constant cross-section, 
whose length is obtained through the flare analysis, the 
only free parameter is the ratio of the diameter of the 
loop's cross section to its total length, or the aspect 
ratio, f3. In the solar context /3 ~ 0.1 is often claimed as 
typical" , altho ugh in the case of act i ve stars higher values 



Maggio et al., 2000| ; [Favata et al., 2000aD 



(up to /3 ~ 0.3, 
have been inferred. 

For a semicircular loop, V = 2Tr For the YLW 15 

flare we thus derive V — 6.1 x 10'^'^-5.5 x 10^'' cm^ 
(where the first value is for /3 — 0.1 and the second for 
(3 — 0.3, as in the following). The density of the flaring 
plasma at the peak of the event can then be derived as 
He = EM /{0.8V), where the factor 0.8 is the scaling 
factor between the number of ions and the number of elec- 
trons for a solar H/He ratio. This assumes a single loop 
uniformly filled with plasma. If the loop is signif- 
icantly filamented (i.e. has a filling factor smaller 
than 1) the resulting density would be a lower limit 
to the actual value. In the prese nt case (using the p eak 
EM = 7.2x10^'' cm-3 reported by flsuboi et al., 2000| ), we 



derive Ue = 3.8 x 10^°-1.3 x 10^" cm~^. The plasma pres- 
sure can then be estimated through = 2nf.kTp = 1600- 
500 dyne cm^^. The equilibrium pressure estimated from 
the peak temperature and the derived length of the flar- 
ing loop through the scaling laws of Rosner ct al. (1978 ) is 
Prtv — 2700 dyne cm~^, showing that, in this case, either 
the actual loop is quite thin, with /3 < 0.1, or that the 
loop is signiflcantly out of equilibrium, perhaps still be- 
ing filled with evaporating plasma when the decay starts. 
The magnetic field necessary to confine the flaring plasma 
is Bq — y/Sirpe — 200-120 G. kG photospheric magnetic 



^ Here, as in the following, the notation x[xl — xu] indicates 
the best-fit value of the quantity x together with the lower and 
upper range of the la confidence interval, which is typically 
not symmetric. 
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fields have been observed through Zeeman splitting in flare 
stars by e.g. Johns-KruU & Valenti (1996). 



3.2. The V773 Tau flare 

V773 Tau is a binary weak-line T Tau star located in the 
Barnard 209 cloud (in the Tau-Aur star- forming region), 
at a distance d ~ 150 p c. The two co mponents have spec- 
tral types K2 and K5 dWelty, 1995| ), with an estimated 
r adius for either com ponent R* ~ 4 Rq = 3 x 10^^ cm 
(Bouvicr ct al., 1995). The interbinary separation is a ~ 
5 X 10^^ cm ( IWclty, 1995| ). 

A large X-ray flare detected during an ASCA obser- 



vation of Barnard 209 has been studied by Tsuboi et al 
1998| ). The decay time-scale of the flare is tlc — 8.2 ks 



with a mELximum observed temperature in the ASCA- 
SIS spectrum of Td ~ 110 MK (i.e. a very ho t event ). The 
peak X-ray luminosity is Lx — 10'^^ 




erg s 



Tsuboi et al. 



(1998 ) note that some of the conditions for the applica- 



bility of a quasi-static analysis are actually violated: in 
particular, the quantity T^ '^^ / EM , which is supposed to 
be constant during the decay, increases toward the end 
of the event, something which [Tsuboi ct al. (19*98 ) inter- 
pret as evidence for reheating. They then use only the 
first part of the decay (in which T^''^^/EM is approx- 
imately constant) to derive the flare parameters. Their 
quasi-static analysis results in a length of the flaring loop 
L ~ 4x 10^^ cm ~ 1.3 i?* and a density n ~ 3x 10^-^ cm"-^. 
The loop is thus derived to be larger than the stars but 
much smaller than the interbinary separation. 

Similarly to the other events discussed here, we have 
analyzed the event using the hydrodynamic m odeling. 



starting from the the spectral analysis results of Tsuboi 



Fig. 2. The evolution of the V773 Tau flare. The top panel 
shows the flare's emission measure (used as a proxy to the 
light curve, e.g. to the X-ray count rate), together with the 
best-fit exponential decay, while the bottom panel shows 
the evolution of the fiare's decay in the logT vs. logn 
plane. The solid line is the best fit to the decay, while the 
dotted line is the path the decay would have followed if 
it had been quasi-static with no sustained heating during 
the decay phase. 



radius i?* ~ 2 Rq = 1.4 x 10^^ cm. A large X-ray flare has 
been detected in its light curve during a ROSAT PSPC ob- 
servation of IC 348, which has been analyzed by Preibisch 



et al. (t998|) . The rCsultej are ejhbwn in Fig, p]: the decay of 



the temperature during the flare is quite shallow, so that 
the shape of the light curve is fully dominated by the tem- 
poral evolution of the heating, with ( = 0.36 ± 0.12 and 
-^(C) = Thc/Tth — 6.4[5.3 — 7.8]. The actual peak tempera- 
ture of the event is very high, at Tmax — 200 MK, and the 
loop length is L ~ 7.5[6. 1-9.0] x 10^" cm ~ 0.25±0.05 i?*. 
The event, notwithstanding its high temperature and X- 
ray luminosity, is originating in a rather small loop, and 
therefore close to the stellar surface. 

Using the same approach as discussed for YLW 15, the 
loop volume can be estimated atV = 2.5 x 10'^^-2.3 x 10'^^ 
cm"^, and the density and pressure can be derived (using 
the peak EM = 5.2 x 10^^ cm-^) at = 1.6 x lO^^. 
5.3 X 10" cm-3 and pe = 2NekTp = 8.6 x 10'*-2.9 x 10* 
dyne cm^^. The equilibrium pressure is prtv = 3.8 x 10"* 
dyne cm~^, compatible with a loop with f3 ~ 0.3. The 
confining magnetic field is estimated at Bq = 1500-850 
G. 

3.3. The LkHa92 flare 

LkHa 92 is a classical T Tau star in the IC 348 star- 
forming complex, with a K spectral type, and an estimated 



et al. (1993|) using the qu asi-static approach. The analysis 
of Preibisch et al. (1993| ) resulted in a decay time-scale 



of the flare tlc = 7.8 ks, and a observed peak temper- 
ature determined from the PSPC spectra Tp ~ 43 MK. 
The peak X-ray luminosity is ~ 5 x 10^^ erg s'-^. The 
quasi-static analysis results in a length of the flaring loop 
L ~ 8 X 10^° cm ~ 0.6 i?, (i.e. a moderate-size loop) and 
a density n ~ 1.5 x 10^"'^ cm^^. 

We have analyzed the event using the approach dis- 
cussed in Reale & Micela (1998), based on a principal 



component analysis, allowing to make optimal use of the 
low-resolution, limited bandpass spectra produced by the 
ROSAT PSPC. The key feature of this approach is the use 
of a non-parametric temperature indicator (the Spectral 
Shape Index, or SSI), rather than a parametric fit to the 
spectrum. 

For the analysis, the flare decay has been binned into 
three intervals, as shown in Fig. |3[ extracting background- 
subtracted spectra for each of them. The quiescent emis- 
sion (determined in the pre-flare phase) was also sub- 
tracted from the flare spectra. The observed decay time 
scale for the event is tlc = 5. 6 [4. 8 — 6.8] ks. The measured 
SSI values have been converted to the intrinsic (unab- 



sorbed) SSI values (as discussed in [Realc fc Micela, 1998]) 
assuming an absorbing column density N{H) = 10^^ 
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Fig. 3. The evolution of the LkHa 92 flare. The top panel 
shows the flare's background-subtracted light curve to- 
gether with the best-flt exponential decay, while the bot- 
tom panel shows the evolution of the flare's decay in the 
SSI (Spectral Shape Index) vs. logn plane (using the 
square root of the count rate as a proxy to the density 
n). The solid line is the best flt to the decay, while the 
dotted line is the path the decay would have followed if 
it had been quasi-static with no sustained heating during 
the decay phase. 



cm^^ (as derived by Preibisch et al., 1993| ). The slope 
in the count-rate vs. SSI plot is = 5.9 ± 1.2. The peak 
SSI corresponds to a maximum temperature for the event 
2max = 110 MK, and the corresponding loop length is 
L = (3.6 ± 1.1) X 10^° cm = (0.26 ± 0.08)i?,, where the 
relatively uncertainty is estimated according to the value 
reported in Table 2 of |Reale fc Micela (199^ ). This flare is 
therefore conflned to a low-lying loop, close to the stellar 
surface, and with a size falling within the range observed 
for solar flares . The heating time s cale th is estimated 
(from Fig. 3 of |Reale fc Micela, 1998|) to be th ~ 2 x tlc- 
The observed flare decay is again dominated by the sus- 
tained heating. Note that the peak temperature implied 
by the peak SSI value is, at llOMK, signiflcantly higher 
than the maximum temperature measured through spec- 
tral fltting (43 MK, see above). This is due to the relatively 
soft bandpass of the PSPC detector. 

Along the same lines discussed above for the YLW15 
event the loop volume can be estimated a,tV = 2.9 x 10"^°- 
2.6 X 10'^^ cm'^ (for (3 — 0.1 and (3 — 0.3 respectively), and 
the density and pressure can be derived (using the peak 



EM 



4.8 X 10^54 cm-3) at He = 1.4 x 10^2-4.8 x 10 
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cm-3 and Pe = 2NekTp = 4.2 x 10^-1.4 x ID"* dyne cm-^. 
The equilibrium pressure is prtv = 1.2 x 10"* dyne cm~^, 
compatible with a /? = 0.3 loop near equilibrium (i.e. flUed 
with plasma). The conflning magnetic fleld is estimated at 
Bo = 1000-600 G. 



Fig. 4. The evolution of the HD 283572 flare. The top 
panel shows the flare's background-subtracted light curve 
together with the best-flt exponential decay, while the bot- 
tom panel shows the evolution of the flare's decay in the 
SSI vs. logn plane (using the square root of the count rate 
as a proxy to the density n). The solid line is the best flt 
to the decay. 



3.4. The HD 283572 flare 

HD 283572 is a relatively nearby WTTS in the Tau-Aur 
star-forming region, at a (Hipparcos-determined) distance 
of 128 pc. Its optical characteristics are described in detail 



by Walter et al. (1987), who derived an MK spectral type 
of G5IV, and a radius of 3.3 i?0, at an assumed distance 
of 160 pc, which becomes 2.7 Rq = 1.9 x lO^^cm at the 
Hipparcos distance. Its X-ray character istics have been 
studied in detail by Favata et al. (1998 ). Thanks to its 
high quiescent X-ray luminosity (Lx — 2 x 10^^ ergs~^) 
it is a well studied object and it has been observed with 
most X-ray telescopes to date. 

A large flare seen with the ROSAT PSPC is rep orted 
(but not discussed in detail) by Stelzer et al. (2GG0| ). We 
have analyzed the event using the same approach based on 
a principal component analysis as used above for LkHa 92; 
note that none of the other flares on PMS flares reported 
by Stclzcr ct al. (2GO0| ) has a sufflcient number of events 
to allow a similar analysis. 

The flare decay has been binned into four intervals, as 
shown in Fig. |^, extracting background-subtracted spectra 
for each of them and subtracting the quiescent emission 
from the flare spectra. Only the flrst three intervals have 
been used for the analysis, given the very low statistics 
of the fourth interval. The observed decay time scale for 
the event is tlc = 12. [9.4 — 18.] ks; the measured SSI 
values have been converted to the intrinsic (unabsorbed) 
SSI values assuming an absorbing column density N{H) 



10 
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(as derived e.g. by Favata et al., 199^ ). The 
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Table 2. The plasma physical parameters derived for the 
flaring events studied here. For each event, the density ng, 
the pressure p and the magnetic field B needed to confine 
the plasma are tabulated under two different assumptions 
for the loop's aspect (length to diameter) ratio, i.e. (3 ~ 0.1 
and P = 0.3. Also, the "equilibrium" pressure prtv that 
the plasma would have in a loop satisfying the scaling laws 
of Rosner et al. (1978) is reported. 
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slope in the count-rate vs. bSl plot ot (,' = 5.1 ± U.S. i'lie 
peak SSI corresponds to a maximum temperature for the 
event Tmax = 48 MK, with a corresponding loop length of 
L= (5.0±1.5)xl0i'^ cm= (0.26±0.08)i?*. The relatively 
uncertainty is also in this case is estimated according t o 
the value reported in Table 2 of |Reale fc Micela (1998|) . 
This flare is therefore also confined in a low-lying loop, 
close to the stellar surface. The heating time sca le th is 
estimated (from Fig. 3 of Reale fc Micela, 1998| ) to be 
TH ^ 2 X TLc ^ 23.4 ks. 

The loop volume can be estimated at y = 7.8 x 10"^°- 
70 X 10^1 cm^ (for ^ = 0.1 and /? = 0.3 respectively), 
and the density and pressure can be derived (using the 
peak EM = 4.3 x 10^'* cm'^) at = 8.3 x 10"- 
2.8 X 10" cm-3 and Pe = 2NekTp = 1.1 x 10'*-3.7 x 10^ 
dyne cm~^. The equilibrium pressure is significantly lower, 
Prtv = 0.8 x 10^ dyne cm~^, thus implying either a loop 
with a very high aspect ratio (3 > 0.3 (which would how- 
ever appear unlikely), or to the Spectral Shape Index un- 
derestimating the peak temperature of the event: the un- 
certainty on the SSI derived for a spectrum of about 1000 
counts (as it is the case for the hotter spectrum in the 
HD 283572 event) is, again from Table 2 of [Reale fc Micela 
(1998|)] Assi — 0.5, so that the la interval for the event's 
peak temperature is 30 - 75 MK, resulting in a la interval 
for the equihbrium pressure Prtv — [0.23.2] x lO'^ dyne 
cm~^, so that the upper value of the la interval is com- 
patible with the P = 0.3 configuration. Also, as shown by 
the other events analyzed here, large flares on PMS stars 
tend to be quite hot, so that indeed the nominal peak tem- 
perature of 48 MK is lower than expected. The confining 
magnetic field is estimated at Bq = 530-300 G. 



4. Discussion 

All flares analyzed in the present paper show evidence for 
strong sustained heating during the flare decay, so that 



the decay light curve is dominated by the temporal heat- 
ing profile and not by the free decay of the heated loop. As 
a consequence, the decay light curve has limited diagnostic 
value for the derivation of the physical parameters of the 
flaring loop, and the quasi-static analysis which has been 
previously applied to the events discussed here invariably 
results in loop sizes which are significantly over-estimated. 
In the presence of such vigorous sustained heating during 
the decay phase, even the sizes derived by hydrodynamic 
modeling (which rely on a parameterization of the heat- 
ing temporal dependence as a decaying exponential) are 
likely to result in upper limits to the act ual loop sizes (as 
demonstrated by the Algol SAX flare, see [Favata, 2000| for 
a discussion). Therefore, the large loop sizes obtained by 
quasi-static analysis (L <; 1 i?*) are in reality likely to be 
significantly smaller, typically by factors ^ 3. 

This result is in line with the results obtained by apply- 
ing the same hydrodynamic modeling approach to other 



classes of active stars, i.e. on Algol (Favata fc Schmitt 



1999; Favata et al., 2000b), on flare stars (Reale fc Micela 



1998; Favata et al., 2000a; Favata et al., 2000c), on the 



zero-age main-sequence star AB Dor (Maggio et al., 2000) 
and on active binaries ( Favata, 2000| ). Sustained heat- 
ing appears to be consistently present during large X-ray 
flares, so that the flaring regions are much smaller in size 
than it would be inferred by assuming free decay of their 
light curve. Therefore the (flaring) corona in very active 
stars appears to be rather low-lying, close to the photo- 
sphere, with no evidence for the very large, extended loops 
which have been deduced in the past from the quasi-static 
analysis of large flares. Pre-main sequence stars are no ex- 
ception to this, as also all flaring events analyzed here 
appear to be due to low-lying structures . 

As discussed by Favata et al. (2000b ), the location of 
large flaring events appears in general to be compatible 
with their being located at high stellar latitude, close to 
the stellar poles, with no evidence for low-latitude, near- 
equatorial structures (as for example typically present on 
the Sun). No large flare has been observed to self-eclipse 
even when its duration is comparable to the stellar rota- 
tional period; given the small size of the flaring regions, 
the polar region is thus the only plausible location for the 
flaring plasma. For the events discussed here, the ones on 
LkHa92, V773 Tau and HD 283572 have durations signif- 
icantly shorter than the plausible rotational period of the 
stars and therefore this line of reasoning cannot be used 
to argue for a polar location of the flaring region. For the 
event on YLW 15, [Montmerle et al. (20001 ) argue that the 
object is a fast rotator, with a period Prot — 16-20 hr. All 
three events reported by Tsuboi et al. (2000| ) have a dura- 
tion of ~ 20 hr each, and in no case there is evidence for 
self-eclipse. Given the relatively small size derived here for 
the flaring region, a polar location for the flaring region 
on YLW 15 appears strongly favored. 

This scenario (a low-lying flaring loop, likely located 
on the polar region of the star) is not compatible with 
the explanation for the lar ge flares on YLW 15 proposed 
by Montmerle et al. (2000[), who, on the basis of the long 
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Table 1. A comparison of the physical parameters derived for the flaring events on PMS stars analyzed here. Td is 
the maximum temperature measured during the flare (indicated only for the events analyzed through a direct spectral 
fitting, and absent for the events analyzed with the Spectral Shape Index - SSI - analysis), and EM the maximum 
emission measure, while Lx is the peak X-ray luminosity, tlc is the decay time for the flare, while L'^^ is the length 
of the flaring loop derived by quasi-static analysis, both in cm and in units of the stellar radius. F(Q is the ratio 
Tlc /Tth between the observed decay time and the natural thermodynamic decay time of the loop (giving a measure of 
the importance of sustained heating in determining the light curve of the decay). Tp is the actual peak temperature in 
the flare derived from the observed temperature Td and the instrument's response or from the value of the SSI, while 
is the length of the flaring loop derived through the hydrodynamic modeling approach (again in cm and in units 
of the stellar radius). Numerical subscripts indicate the power of 10 by which the relevant quantity has been scaled. 

Star Instr. Td EM54 tlc Lf^ ngf F(C) Tp Lfi^ 

MK cm"'^ ks cm {R,) cm'^ MK cm (fi«) 

YLW15 ASCA-GIS 64 15 37^9 100 (3.5) 5 2?r 150 46.0 (1.5) 

LkHa92 ROSAT-PSPC 43 50 5.6 8 (0.6) 15 4.3 110 3.6 (0.3) 

V773 Tau ASCA-GIS 110 100 8.2 40 (1.3) 30 6.4 200 7.5 (0.3) 

HD 283572 ROSAT-PSPC - 4 12.0 - - 4.6 48 5.0 (0.3) 



loop size derived from the quasi-static analysis argue for 
a flaring loop with one foot anchored on the stellar sur- 
face and the other onto the (edge of the) accretion disk. 
The smaller loop size for the YLW 15 flare implied by the 
hydrodynamic modeling discussed in the present paper 
makes such loops implausible, and therefore it also casts a 



doubt on the framework used by Montmerle et al. (2000) 



to explain the repeated flaring events seen on YLW 15 (i.e. 
the stressing of magnetic field lines spanning both the star 
and the disc). 

A summary of the characteristics of the event for each 
of the flares studied here is reported in Table |l], and a 
summary of the resulting plasma parameters in Table ^j. 
For the events on the CCTS LkHa 92 and on the WTTS 
V773 Tau and HD 283572 the resulting loop sizes are con- 
sistently around L ~ 0.3 i?,. The event on the proto- 
star YLW 15 is somewhat larger, with a nominal length 
L ~ 1.5 However, the pressure derived even assuming 
/3 = 0.1 is significantly lower than the equilibrium pres- 
sure derived from the [Rosner et al. (1978 ) scaling laws. 



This points toward the loop being far from equilibrium 
and only partially filled with hot plasma, so that the pro- 
cedure used here is likely to over-estimate the actual loop 
size. Both these elements make us consider the length es- 
timate derived here as an upper limit to the actual size 
of the fiaring region. Therefore the characteristics of the 
X-ray fiare observed on YLW 15 appear to be similar to 
the ones observed on older stars, so that there may not 
be any need to invocate a separate mechanism to explain 
the observed behavior. Whether flares on protostars are 
in general more "out of equilibrium" than flares on older 
stars is of course impossible to say on the basis of a single 
event, and clearly study of a significant number of flaring 
events will be needed to assess whether any such general 
difference exists. 
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